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MODELING THE SURFACE X-RAY EMISSION AND VIEWING GEOMETRY OF PSR J082 1-4300 IN PUPPIS A 
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ABSTRACT 

We present a model for the unusual X-ray pulse profiles of PSR J082 1-4300, the compact central object in 
supernova remnant Puppis A. We show that a pair of thermal, antipodal hot-spots on the neutron star surface is 
able to fully account for the pulsar's double blackbody spectrum and energy-dependent pulse profile, including 
the observed 180° phase reversal at » 1 .2 keV. By comparing the observed pulse modulation and phase to the 
model predictions, we strongly constrain the hot-spot pole (£) and the line-of-sight (ip) angles with respect to 
the spin axis. For a nominal radius of R = 12 km and distance D = 2.2 kpc, we find (£,tp) = (86°, 6°), with 
1(7 error ellipse of (2°, 1°); this solution is degenerate in the two angles. The best-fit spectral model for this 
geometry requires that the temperatures of the two emission spots differ by a factor of 2 and their areas by a 
factor of ~ 20. Including a cosine-beamed pattern for the emitted intensity modifies the result, decreasing the 
angles to (84°, 3°); however this model is not statistically distinguishable from the isotropic emission case. We 
also present a new upper limit on the period derivative of P < 3.5 x 10~ 16 (2a), which limits the global dipole 
magnetic field to B s < 2.0 x 10 11 G, confirming PSR J0821-4300 as an "anti-magnetar." We discuss the results 
in the context of observations and theories of nonuniform surface temperature on isolated NSs of both weak 
and strong magnetic field. To explain the nonuniform temperature of PSR J082 1-4300 may require a crustal 
field that is much stronger than the external, global dipole field. 

Subject headings: pulsars: individual (PSR J082 1-4300) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

The young X-ray pulsar PSR J082 1-4300, associated with 
the Puppis A supernova remnant (SNR), is one of three 
pulsars i n SNRs that are spinnin g down nearly imper- 
ceptibly dGotthelf & Halpernl 120091 Paper I). The age of 
PSR J082 1-4300 is 3.7 ky r based on the prop er motion of 
oxygen knots in Puppis A (Winkl er et alj|1988b . and its dis- 
tance is 2.2 kpc fro m H I absorption features to the SNR 
dRevnoso et al.|[l995l) . In the context of the magnetic dipole 
model, these pulsars are a new physical manifestation of neu- 
tron stars (NSs), "anti-magnetars" born with weak magnetic 
fields possibly related to slow natal spin (Gotthelf & Halpern 
2008). They were drawn from the previously defined class of 
central compact objects (CCOs) in SNRs, which are charac- 
terized by their steady, predominantly thermal X-ray emis- 
sion, lack of optical or radio counte rparts, and absenc e of 
pulsar wind ne bulae (see reviews by iPavlov et al J 120041 and 
iDeLucal 120 08). Currently, 7-10 objects are known or pro - 
posed to be CCOs (for a list, see lHalpern & Gotthelfll2010h . 
and are therefore candidates for anti-magnetars. 

The XMM-Newton discovery observations of 
PSR J082 1-4300 revealed an abrupt 180° phase rever- 
sal of its quasi-sinusoidal pulse profile at an energy of around 
1.2 keV (Paper I). The X-ray spectrum of PSR J082 1-4300 
was fitted with a two-temperature blackbody model, both 
temperatures being seen at all rotation phases, while the 
cross-over energy of the spectral components coincides with 
the energy where the pulse reverses phase. These detailed 
properties afford the opportunity to construct a highly 
constrained model of the NS surface emission geometry. 

In this Paper, we present a quantitative verification of the 
geometrical model for PSR J082 1-4300 proposed in Paper I. 
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By reproducing the detailed pulse profile behavior, we are 
able to specify the surface emission areas and viewing geome- 
try to within < 2°. Our treatment includes general relativistic 
effects of light deflection and gravitational redshift, and exam- 
ines the effects of local anisotropy (beaming) in the emitted 
radiation. We describe the antipodal hot-spot model in Sec- 
tion 2, compare the energy-dependent modulation to the data 
for the range of allowed geometries in Section 3, and explore 
if it is possible to limit the neutron star radius in conjunction 
with the estimated distance. In Section 4 we discuss some im- 
plications of the model results, and in Section 5 compare with 
pulsars of other types. 

2. THE EMISSION MODEL 

Our method of modeling the emission from spots on 
the surface of a N S follows the derivation given by 
iPechenick et al.l(ll983l). w ith some generalizations introduced 
bv lPerna & Gotthelf! (120081) . The radiation comes from a hot 
spot of blackbody temperature 7), and angular radius /3k, and 
an antipodal warm spot of lower temperature T w and angular 
radius /3 W . The remainder of the surface is assumed to be at a 
uniform temperature 7ns < T w . The geometry is indicated in 
Figure Q] 

We use 7(f) to indicate the phase of rotation instead of the 
common notation 0(f)- Phase 7 = corresponds to the clos- 
est approach of the hot spot to the observer, while the phase 
of rotation is related to the angular rotation rate of the star 
through 7(f) = fit . We indicate with otySf) the angle that the 
hot-spot axis makes with the line-of-sight. a/,(0 is a func- 
tion of the angle £ between the hot-spot axis and the rotation 
axis and the angle ip between the line-of-sight and the rotation 
axis, by means of the relation 

cth(t) = arccos[cos-0cos£; + sin?/' sin £ cos 7(f)] . (1) 

For each set of angles £ and ip, the angle a w that the axis of the 
opposing warm spot makes with respect to the line-of-sight is 
simply a w (t) = ir-ah(t). 
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The spherical coordinate system (6,<p) is defined with re- 
spect to the line-of-sight as the z-axis. Due to general rela- 
tivistic effects, a photon emitted at a colatitude 9 reaches the 
observer only if emitted at an angle 5 with respect to the per- 
pendicular to the NS sur face. The two angles are related b y 
the ray-tracing function 3 (Pec henick et al.l ll983: Page 1995) 
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0(6)= / xdu 
Jo 

(2) 

having defined x = sin 6. Here, R/R s is the ratio of the NS 
radius to Schwarzschild radius, R s = 2GM/c 2 (we will assume 
M=\AM Q ). 

The hot spot is bounded by the conditions: 
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On the other hand, it is identified through the condition 

9 <9 l Ka h ,/3h,4>) if and f3 h > a h , 



(6) 

where the outer boundary f^(a/,, /?/,, </>) of the spot is computed 
by numerical solution of the equation 

cos/3/, = sin#J sina/, cos</>+cos6^ cos a/, . (7) 

The antipodal warm spot is described on the surface of the 
star through the same conditions, but with the substitutions 
Ph — > Ay and a h — > a w . 

We assume that the emission from the hot and warm spots 
is blackbody, of uniform temperatures 7), and T w , respec- 
tively. The spectral function is then given by n(E,T) = 
l/[exp(E /kT)- 1], where the temperature T{9,<fi) is equal 
to 7), or T„ if 9 and are inside either of the spots, respec- 
tively, and it is equal to zero outside. Given that the presence 
of NS atmospheres and their elemental composition are yet 
to be firmly established, here we first model isotropic radia- 
tion. Then, we explore how the obtained constraints depend 
on the assumption of forward beaming of the radiation, as 
found in magnetiz ed, light element atmosphere models (e.g., 
iPavlov et al.lll994l) . by approximating beaming 4 as intensity 
1(8) oc cos 5 n . The observed spectrum as a function of phase 
angle 7 is obtained by the standard method of integrating the 
local emission over the observable surface of the star, ac- 
counting for the gravitational redshift of the radiation follow- 
ing HagidHH: 

x / ^k ) (9,4>)n[E OQ e-^-T(9,4 > )] (8) 
Jo 27r 

3 To improve the comp utational efficiency of this equation we use the ap- 
proximation presented in Beloborodov (2002). 

4 This approximation is based on the assumption that the magnetic field 
is normal to the surface in the heated regions; this is plausible since heat 
transport along the B field lines is enhanced in the outer envelope, at least for 
B > 10 10 G. See §4 for further discussion of this issue. 
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FIG. 1. — Emission geometry on the surface of the NS for the model pre- 
sented herein: a hot spot of temperature 7), and angular size ft, and an an- 
tipodal spot of temperature T w and angular size /9 W . As the NS rotates with 
angular velocity f2, the angle a(t) is a function of the phase angle 7(f) = Qt, 
the angle £ between spin axis and hot spot axis, and the angle tji between spin 
axis and line-of-sight. 



in units of photons cm" 2 s" 1 keV -1 . In equation ([8]), the NS 
radius and photon energy as observed at infinity are given by 
= Re~ A ' and Eoo = Ee A ', where E is the energy emitted at 
R, and A s is defined as 



R 



(9) 



The phase-averaged flux is computed as F a vg(£oo) = 

In addition to the basic two-temperature antipodal spot 
model described above, the spectrum of PSR J082 1-4300 re- 
quires an additional narrow line-like component around 0.77 
keV, possibly an electron cyclotron feature in emission (Pa- 
per I). Furthermore, as shown in Paper I, this spectral feature 
is associated exclusively with the larger spot, of temperature 
T w . In the current study, we include the best fitted Gaussian 
line as an additive component to our basic model. With no 
other information about its spatial distribution, this emission 
is assumed to be spread uniformly over the surface of the 
warm spot only. This line emission is shown to account for 
a notable increase in the observed modulation below 1 keV, 
as described in Section 3. 

3. MODELING THE ENERGY-DEPENDENT MODULATION 

The surface emission geometry of PSR J0821-4300 is 
highly constrained by its unique energy dependent pulse pro- 
file. As shown in Paper I, the quasi-sinusoidal signal has a 
background subtracted pulsed fraction of « 11% in the en- 
ergy band 0.5-4.5 keV, with an abrupt 180° change in phase 
at 1.2 keV, around which the modulation evidently cancels 
out. This behavior is indicative of a geometry having the 
symmetry of Figure [2 namely, a pair of antipodal spots of 
different temperatures. Our goal is to match the observed 
pulse profile (modulation and phase) in three interesting en- 
ergy bands, 0.5-1.0, 1.0-1.5, and 1.5-4.5 keV, using the an- 
tipodal model, by exploring the range of all possible viewing 
and hot-spot geometry pairs (£,,ip; see Figure [T]), and fitting 
for the correct one. 
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We summarize our proceedure as follows. We started by fit- 
ting the X-ray spectrum of Paper I using the antipodal model 
to compute the temperatures and spot-sizes corresponding to 
all spot and viewing angles £ and ip. We then used these mod- 
els to compute predicted pulse profiles in the three bands as 
a function of (£,ip), to compare with the observed profiles. 
Clearly, only certain geometries will produce a phase shift; 
for example, there will be no shift if the spot axis is nearly 
coaligned with both the viewing direction and spin axis, as 
only one spot remains in view as the star rotates. Finally, we 
considered the effect of radiative beaming, and repeated our 
analysis for a range of NS radii. In the following we describe 
our procedure and results in detail. 

For a given set of viewing angles the antipodal flux model 
can be integrated over phase to provide a direct compari- 
son with the observed spectra. We have incorporated equa- 
tion ^ into an "add itive model" fo r use in the XSPEC spec- 
tral fitting software dArnaudlfl996l) . The coded model com- 
prises 13 parameters: the NS radius and distance (R,D), three 
blackbody temperatures (7, v , 7/,, 7ns), two spot angular sizes 
Ww, Ph), two geometrical angles (£, ijj), the rotation phase (7), 
and the Gaussian emission-line center, width, and flux. The 
column density is fixed at the best value determined in Paper I 
from a fit to the overall spectrum, Nh = 4.8 x 10 21 cm" 2 . In the 
spectral fits, the normalization is set to unity so that the flux is 
determined by R and D, and implicitly takes into account all 
relativistic effects noted in Section 2. 

The model allows for a uniform temperature 7ns in the 
inter-spot area, but since none is necessary for an accept- 
able fit to the spectrum or pulse profiles, we set 7ns = 0. 
Nevertheless, to place a model-independent upper limit on 
7ns, we have simulated blackbody spectra in XSPEC for 
R = 12 km and D = 2.2 kpc (Reynoso et al. 1995), increas- 
ing the temperature until the model exceeds the spectrum of 
PSR J0821-4300 at low energies. Since the model counts 
depend on the interstellar col umn density, we ass umed here 
the largest value measured by iKatsuda et alJ (1201 Ol) for fila- 
ments in the Puppis A remnant in Chandra data, Nh = 5.5 x 
10 21 cm" 2 . This yields a conserv ative 3er upper limit of 
7^ < 0.15 keV; iHwang et all (120081) report a significantly 
lower value of A^h = 3 x 10 21 cm" 2 from Suzaku measure- 
ments, which would allow a smaller limit on T-gg, if appli- 
cable. 

To map the antipodal model for PSR J082 1-4300 as a func- 
tion of (£, t/0 we generated best fit model parameters over the 
grid spanning (0° < £ < 90°; 0° < if) < 90°) in steps of 10° by 
systematically fitting for T w , Th, f3 w , and /?/,. Since we fitted 
the phase-averaged spectrum (presented in Paper I and de- 
scribed therein), a reasonable approximation for the equiva- 
lent phase-averaged model flux is F avg = [FQy = 0°) + F(7 = 
180°)]/2. The flux of each blackbody component depends on 
the size of its respective emission spot as parameterized by 
its (3, which sets the normalization for that component. As 
we consider the range of viewing orientations, the size of the 
spots must be adjusted to keep the model flux fixed to that of 
the observed value. This is because the projected flux from 
the phase integrated emission strongly depends on the view- 
ing geometry (£,V0> which is not known a priori. In compar- 
ing with data, for each trial (£, ip) pair, we fit for the spot sizes 
(PwiPh) that correctly normalize the blackbody components 
to match the observed spectrum. 

The grid of best fitted models with values of (f3 w , /?/,, 7, v , 7),) 
describes the flux from PSR J082 1-4300 as a function of en- 



ergy and phase (FlEoo,^]) according to equation <[8J f° r an y 
realizable geometry 5 , the latter parameterized by (£,V>) and a 
set of fixed R,D, and emission-line parameters. We find that 
the spot temperatures remain almost constant in (£, ip) and use 
the average value given in Table 1 in our final fits. The result- 
ing fitted values of /3 W and /?;, for the case R = 12 km are inter- 
polated to 1 x 1 pixels on a 90° x 90° grid and displayed in 
Figure |2] 

We are now prepared to compute the set of phase-resolved 
modulations in the three interesting energy bands, which 
manifest the phase shift. Specifically, the pulse in the soft 
0.5 - 1 .0 keV band peaks at rotation phase 7 = 180°, the hard 
1.5-4.5 keV band peaks at 7 = 0°, while the modulation 
in the intermediate 1.0-1.5 keV band cancels out to < 1%. 
The comparison is best done in counts space rather than flux, 
to make use of the counting statistics for errors. The back- 
ground counts are measured for each band and added to the 
model counts. For each (£,ip) pair, the antipodal flux model, 
including interstellar absorption, is folded through the XMM- 
Newton EPIC response function used for the spectral fits (see 
Paper I). For each energy band, we record the magnitude of 
the model modulation f p defined as 

F^-NF^ ^ (1Q) 

Aot 

where the flux F tot is the total flux in the band, N is the num- 
ber of bins in the light curve (ten in this case), and N F m i n is 
the "unpulsed" flux in the band, determined from the bin of 
minimum flux. The results are shown in Figure[3]for the soft 
and hard energy bands. 

For a given (£, iji) pair, the model profile in ten phase bins is 
compared to the observed counts in the three bands, for a total 
of 30 bins, using the \ 2 statistic. Since the total count nor- 
malizations in each band were fixed earlier by the spectral fits 
via the /3 parameters, the mean counts of the models match 
the means of the three observed light curves. This reduces 
the number of degrees of freedom by three, from 30 to 27. 
The (£, rp) region that best matches the measured pulsed frac- 
tion in the soft and hard bands, f p = 11% ± 2% corrected for 
background, is indicated by the contours in Figure [3] While 
there is a range of allowed angles in each band separately, they 
overlap in only two small regions corresponding to the mini- 
mum x 2 . This uniquely identifies the geometry of the system, 
within statistical uncertainty. More specifically, we find that, 
for R = 12 km, the most likely geometry is specified by the 
combinations (ip,0 = (86°, 6°) or (6°, 86°). Since the flux de- 
pends on the angles £ and ip only through the parameter a/, in 
equation (JTJ, it is symmetric with respect to an exchange of 
£ and yielding the two possible solutions. This is evident 
in the \ 2 ma P that compares the model and observed profiles 
in the three bands (Figure 0J. In one solution the spin axis 
is nearly perpendicular to the line of sight, while in the other 
solution the two lines are nearly parallel. Figure [5] shows the 
fit to the pulse profiles for these equivalent best solutions. Fi- 
nally, the parameters for the model that resulted in the small- 
est x 2 for eac h test radius R are given in TableQ] In all cases, 
the reduced % 2 is near unity, indicating an excellent match to 
the observed pulse profiles. 

Figure [6] provides a graphic explanation of the origin of 
the observed energy-dependent modulation and phase-shift 

5 Not all values of (£,,ip) can fit the spectrum, as the available flux in some 
cases, e.g., £,ij> both near zero, are insufficient to match the data even with fi w 
at its maximum value of 90° . Such cases produce unacceptable \ 2 statistics. 
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FIG. 2. — Best fitted values for the sizes of the two emitting spots on the surface of PSR J0821-4300 as a function of geometry parameters for a NS 

radius of R = 12 km. Left: Map of warm-spot size parameterized by its angular radius /3„ ... Right: Corresponding map for hot-spot angular radius /3;,. The spot 
size is computed at intervals of 10° in i/i and £, and interpolated to 1°. Note the very different size ranges of each spot. Using these results, the modulation as a 
function of viewing geometry is then computed and compared to the data, constraining the allowed geometry as shown in Figure [3] This procedure is repeated 
for a range of NS radii, each requiring a new pair of maps, with results given in Tabled 
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FIG. 3 . — The model modulation derived in two energy bands as a function of geometry parameters (£, ip), for a NS radius of R = 12 km. The pulsed fraction 
ranges from 0-65% and is scaled linearly, with yellow denoting the largest modulation. The solid line in each panel indicates the contour of the measured pulsed 
fraction in the soft energy band (left panel) and hard energy band (right panel), while the dashed lines give the lcr error range. The possible geometry is then 
strongly constrained by the intersection of the two contour regions, as shown in Figure|4] Note that both antipodal spots contribution to the modulation shown in 
each of these plots. 



seen from PSR J082 1-4300. For the set of model parameters 
that best fit the data, we graph separately the phase-averaged 
fluxes for the two spots to gauge their contribution to the light 
curve. At the lower energies, the large warm spot dominates 
the spectrum and the light curve peaks when this component 
is in view, while at the higher energies, the small hot spot 
dominates and it peaks in view 0.5 cycles later. Thus, the 
energy at which the dominant spectral component switches 
is around 1.3 keV, and the peak phase necessarily shifts by 
180° at this energy due to the antipodal symmetry in longi- 
tude. This phase reversal, at this energy, is thus a direct conse- 
quence of the crossing of the spectral components of different 



temperatures. The agreement between the energy of phase re- 
versal and the energy of the spectral cross-over point provides 
direct evidence of the correctness of the model. 

Figure[6]also shows that the predicted modulation is 19.6% 
at the highest energies, where 100% of the flux comes from 
the Th blackbody component. At lower energies, the modu- 
lation is reduced exactly in proportion to the increased con- 
tribution of the T w or 7), spot, depending on which one domi- 
nates the flux at a given energy, f p (E) = 19.6% x \F(E;T W )- 
F{E; Th)\/[F(E; T W ) + F(E; 7),)]. For our broad energy bands, 
used to compare the model with data, the observed modula- 
tion is evidently weighted by the total flux over the band. This 
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FIG. 4. — Contours of obtained by comparing the pulse profiles of the 
antipodal model with the data in three energy bands, for a range of angles 
£ and ip, as described in the text. The lo\ 2o\ and 3cr confidence levels 
are shown for the best match for a NS radius of R = 12 km. The results are 
degenerate with respect to an interchange of £ and ip. The minimum xi f° r 
the viewing geometry parameters is obtained at ip = 86° and £ = 6° , evidently 
providing a strong constraint. The geometries that manifest a phase shift (or 
not) are indicated. 
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PSR J0821-4300 PHASE-AVERAGED MODEL SPECTRUM 

N H = 4.8x10 z 'cm 
R = 12 km 
D = 2.2 kpc 
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FIG. 5. — Pulse profiles generated using the antipodal model (solid line) that 
best matches the observed data (histogram) across the three selected energy 
bands for PSR J082 1—4300. The parameters of this model assuming a NS 
radius R = 1 2 km are given in Table [T] The measured background in each 
band has been added to the model to enable a direct comparison. 



is especially clear in the highest energy band (1.5-4.5 keV), 
in which most of the photons are from the lower end of the 
band where the modulation is significantly less than 19.6%. 
The modulation in the middle band, which spans the phase 
shift with equal counts, is mostly canceled out. That a similar 



FIG. 6. — The antipodal explanation for the observed energy-dependent 
modulation and phase shift between the soft and hard X-ray bands seen in 
PSR J0821-4300. Top: Contribution of each phase-averaged model spectral 
component to the two-blackbody model, plus a Gaussian emission line. The 
dotted line shows the warm blackbody component without the emission line 
contribution. Bottom: The model pulse modulation as a function of energy. 
The strong energy dependence results from the relative contributions of (out 
of phase) overlapping (in energy) flux from the two spectral components. The 
maximum modulation is 19.6% at the highest energy. However, the modu- 
lation in Band 3 is smaller, as it is weighted by lower-energy photons. The 
emission-line contribution is evident as the difference between the solid and 
dotted lines. The phase of the peak of the light curve follows the dominant 
spectral component at each energy, and is restricted by the symmetry of the 
model to either or 0.5 cycles. Where the spectral components cross, the 
phase must shift by 180° as observed (see Figure|3]and Paper I). 



pulsed fraction is measured in the lower and higher bands is 
largely coincidence in this case. Figure [6] also shows the con- 
tribution of the spectral line at 0.77 keV; a small but signifi- 
cant increase in the modulation results. The photon statistics 
of the current data do not allow a more detailed comparison 
of the modulation as a function of energy. 

Next we explored the dependence of the fitted parameters 
on the radius R of the NS. As discussed above, due to general- 
relativistic effects R is not simply a normalization of the flux; 
this requires us to compute j3 maps for each test value of R. 
Therefore, we repeated our full analysis for the range of val- 
ues 9 < R < 14 km, in 1 km increments. These results are 
presented in Table Q] In principle, the relativistic effects can 
lead to a preferred R, but statistically, no unique radius is 
suggested here. The general trend is that the angles j3 w ^p3h, 
and £ increase with decreasing radius. Two counteracting ef- 
fects, both due to flux conservation, influence the spot sizes. 
Gravitational redshift decreases the inferred emission area on 
a more relativistic (smaller) star to compensate for the larger 
needed surface temperature. On the other hand, for a fixed 
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distance D, the spot angular size increases on the smaller star. 
For the values of the fitted parameters here, the latter effect 
tends to dominate over the former, reducing the modulation 
from a smaller star for the same viewing geometry. How- 
ever, the most important effect influencing the amplitude of 
the modulation is the gravitational deflection of light rays, 
which acts to suppress the pulsed flux for the smaller star. 
Therefore, in order to reproduce the same observed level of 
modulation, a smaller R requires a larger modulation of the 
viewing angle a(t), which, according to equation (fTJ, is ob- 
tained by increasing either £ or ip. This explains the trends in 
Table QJ 

We then repeated our analysis for a locally anisotropic in- 
tensity pattern, I(S) oc cos S. For the nominal radius R = 
12 km, the combined constraints from the pulsed fraction 
and the phase shifts moves the best fitted viewing angles to 
(£, ip) = (84°, 3°), with a larger uncertainty than found for the 
unbeamed case. This behavior is readily accounted for. For- 
ward beaming enhances the emission in the direction of the 
spot axis, hence increasing the differences in observed flux 
as the axis of the spot moves toward and away from the ob- 
server. As a result, a fit to the observed level of modulation 
requires smaller values of the viewing angles. The pulse pro- 
file is found to remain sinusoidal and is statistically indistin- 
guishable from the isotropic intensity case. For the specific 
geometry of PSR J082 1-4300, such beaming is a weak ef- 
fect in narrowing the pulse profile because of our relatively 
unmodulated views of the two antipodal spots at glancing an- 
gles. 

In the antipodal model, the energy dependent phase rever- 
sal is a direct consequence of switch in dominance between 
the two blackbody spectral components of different tempera- 
tures. As mentioned previously, this agreement between the 
cross-over energy in phase and spectra is distinct feature of 
this symmetric model, where the emission spots are exactly 
antipodal. However, similar light curves can be obtained if the 
spots remain opposite in longitude, but are allowed to move 
closer in latitude. Such geometry can still produces a 180° 
phase reversal. Allowing this additional degree of freedom, 
the angles (£,ip) would not be so strongly constrained as in 
the antipodal model. More generally, if the spot locations dif- 
fer in longitude by A7 < 180°, such an "offset" model allows 
the possibility of a continuous phase shift as a function of 
energy. The effect of asymmetric spot locations is clear in 
the energy-dependent model profiles, but any such effect in 
PSR J082 1-4300 is not apparent, and would require higher 
quality data to discern. Description of the geometry of off- 
set models can be found in Bogdanov et al. (2007, 2008), for 
example. 

4. DISCUSSION OF MODEL RESULTS 

Using an antipodal spot model, we have accounted for 
all of the details of the XMM-Newton observations of 
PSR J082 1-4300 described in Paper I. In particular, we can 
reproduce the overall spectral shape, energy-dependent pulsed 
modulation, and abrupt 180° phase reversal at the cross-over 
energy of the fitted blackbody components. In so far as 
no observed phenomena remains unmodeled, and no unob- 
served features are predicted, the antipodal model provides 
a credible description of the geometry of emission from the 

6 The larger of the two angles (£, tp) is unchanged as a function of radius 
relative to its 2° error, while the smaller one is clearly decreasing with radius 
as compared to its 1 error. 



TABLE 1 

Model Results as a Function of NS Radius 



Parameter 


Unc. fl 


9 


10 


R(km) 
11 


12 


13 


14 


kT„ (keV) 


3% 


0.29 


0.28 


0.27 


0.26 


0.26 


0.25 


1cT h (keV) 


3% 


0.57 


0.54 


0.53 


0.51 


0.50 


0.50 


A* 


6% 


39° 


37° 


35° 


34° 


32° 


31° 


0h 


11% 


8° 


7° 


7° 


7° 


6° 


6° 


£ or ip 


2° 


86° 


87° 


87° 


86° 


87° 


87° 


ip or £ 


1° 


90 


7° 


7° 


6° 


5° 


5° 


Xl (27 DoF) 




0.94 


1.00 


0.99 


1.00 


0.95 


0.97 



NOTE. — Spectral fits with fixed parameters Nh = 4.8 X 10 21 cm -2 
and D = 2.2 kpc. Gaussian line model parameters are fixed at en- 
ergy E°° = 0.77 keV, width a = 0.05 keV, and flux normalization 
1.8 x lO" 4 phcnr 2 s _1 . 

" The lcr uncertainties in the spectral parameters are estimated by run- 
ning the XSPEC error command; the uncertainties in £ and 1/1 are de- 
termined from the x 2 raa P of Figure [4] 



CCO in Puppis A. The full data set can be reproduced, 
with slight differences in the best fitted parameters, assum- 
ing either isotropic or forward-beamed emission. Differen- 
tiating between these assumptions will require observations 
with higher statistics. 

By matching the observed modulation in three broad en- 
ergy bands, we are able to restrict the angles that the hot-spot 
axis and the line of sight make with respect to the spin-axis 
to within < 2°, up to the degeneracy between these two an- 
gles. Either the spin axis lies nearly parallel (6°) to the line- 
of-sight, with the hot-spot axis at 86°, or the hot-spot axis is 
nearly co-aligned with the spin-axis, but perpendicular to the 
line-of-sight. In the absence of a strong physical motivation 
to prefer one of these configurations over the other, we note 
that the a priori probability of the spin axis lying 6° ± 1 from 
the line of sight is only 3.6 x 10~ 3 , while the probability that 
it is at 86° ±2° is 7.0 x 10" 2 , a factor of 20 larger, although 
still small. We note that the specific orientations that fit the 
observations of PSR J0821-4300 are not the only ones that 
allow phase reversals in the two-temperature model. Rather, 
phase reversals are found in the majority of configurations of 
Figure |4] 

Comparing with the other CCO pulsa rs, we see that PSR 
J 1852 +0040 in the SNR Kesteven 79 (lHalpern & Gotthelfl 
120101) also has a two-temperature X-ray spectrum, but its 
highly modulated pulse (f p = 64%) is single-peaked and 
virtually invariant with energy. Because of this, its emit- 
ting regions are likely to be concentric, or nearly so. In 
the case of IE 1207.4-5209 in PKS 1209-51/52, there are 
large variat ions in pulse phase a nd amplitude as a fu nction 
of energy (Pa vlov et al.l [2002b; D e Luca et al.1 [2004), with 
the largest pulsed fraction coinciding with the strong ab- 
sorption lines in the unique spectrum of this pulsar. This 
effect may be a manifestation of angle-dependent scat- 
tering in cyclotron lines, which is the favored identifica- 
tion of the spectral features considering the upper limit of 
B s < 3.3 x 10 11 G on the surface dipole field from the ab- 
sence of spin-down (iGotfhelf & Ha lpern 2007). The data on 
IE 1207.4-5209 should be fitted with detailed atmosphere 
models that include quantum t reatment of the cyclotron har- 
monics (ISuleimanov et al.ll2.010h . 

Applying this model to surface thermal emission from 
CCOs with weak magnetic fields (anti-magnetars) is an espe- 
cially apt use, in that additional complicating emission mech- 
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anisms that are evident in other classes of pulsars (see Sec- 
tion 5) are absent in CCOs. Such contributions include non- 
thermal magnetospheric emission in spin-powered pulsars, 
polar-cap heating from backflowing particles, and transient 
and variable heating from magnetic field decay in magnetars. 
The first two extra contributions can be significant even for 
recycled millisecond pulsars, which are now known to be ef- 
ficient 7-ray emitters (Abdo et al. 2009a, b). Observations and 
upper limits on spin-down of CCOs indicate spin-down lu- 
minosities that are smaller than their thermal X-ray luminosi- 
ties, and dipole magnetic fields of order 10 10 " 11 G, remarkably 
small for young pulsars. These properties imply that none 
of the above-mentioned emission and surface heating mech- 
anisms can be significant, and constrain the effects that may 
be responsible for the multiple temperatures that are a ubiqui- 
tous feature of CCO spectra, even those that have not yet been 
observed to pulse. 

For the assumed distance to PSR J082 1-4300 of 2.2 kpc 
and a radius of 12 km, the best match for the modulation 
fixes the extent of the hot and warm regions to angles fth = 
6.°6±0.°5 and /3 W = 34.°0± 2°, representing 0.33% and8.5% 
of the surface area, respectively. The existence of a hot spot is 
difficult to understand in the context of a weakly magnetized 
NS, as it requires a mechanism to confine the heat to such a 
small area. Using XMM-Newton, we obtained a new period 
measurement of PSR J0821-4300 on 2010 May 2 using the 
identical observational setup and analysis as described in Pa- 
per I. We obtained another measurement of the pulsations on 
2010 Aug 16, from a Chandra CC-mode observation. These 
results will be presented in a future publication. The period 
is found to be unchanged from the values observed in 2001. 
In combination with the previous measurements listed in Pa- 
per I, this places a 2er limit on P of < 3.5 x 10~ 16 and, under 
the assumption of dipole spin-down, B s < 2.0 x 10 11 G, con- 
firming PSR J082 1-4300 as an anti-magnetar. Given the cor- 
responding upper limit on spin-down luminosity of < 1 x 10 34 
erg s , the hot-spot luminosity of w 2 x 10 33 erg s" 1 can 
hardly be attributed to external heating by backflowing par- 
ticles. The same problem was discussed in the context of 
the highly pulsed emission from PSR J 1852+0040 in Kes 79 
(lHalpern & GottheldloToh . 

Possible explanations for the properties of CCOs are largely 
focused on magnetic field induced aniso tropies in the surface 
tempe rature of a NS, as proposed by iGreenstein & Harikel 
( 1983), in which strongly enhanced conductivity in the direc- 
tion parallel to the magnetic field is matched by a correspond- 
ing reduction in the perpendicular direction. The effect of the 
magnetic field on the heat transport of the crust and envelope 
of neu tron stars has been invest i gated b y a number of authors 
(e.g. iHeyl & Hernquistl 1 19981 120011: [Potekhin & Yakovlevj 
l200lHLail200lHGeppert et al.l2004ll2006tlPerez-Azorfn et all 
l2006at iPons et al.1 120091) . While heat transport in the core 
(p ~ 1-6 x 10 14 g cm" 3 ) is expected to be roughly isotropic 
due to proton superconductivity, anisotropy of heat transport 
becomes pronounced in the outer envelope (p < 10 10 g cm" 3 ) 
for field strengths B > 10 10 G, and it extends deeper into the 
whole crust for higher fields, B > 10 12 - 10 13 G. The main 
question is whether subsurface fields in CCOs can be strong 
enough to affect heat transport to the extent required, while 
not exceeding the weak external dipole field as constrained 
b y their spin-down prop ertie s . 

Gep pert et al.l (120041) discussed the differing effects of a 
poloidal magnetic field in the core of the NS, versus one 



confined to the crust, the true configuration being a matter 
of uncertainty. From a core field, any surface temperature 
anisotropy is expected to be small, while a tangential crustal 
field insulates the magnetic equator and conducts heat to the 
magnetic poles. A tangential crustal field may be indicated 
for CCOs, because it can lead to small hot regions where 
the field emerges normal to the surface, while contributing 
very littl e to the external dip ole field. Of particular inter- 
est here, Gepp ert et al.1 d2006) found that, if the crustal field 
consists of both a dipolar poloidal and a toroidal compo- 
nent, then configurations can be realized in which two warm 
regions of different sizes are separated by a cold equatorial 
belt. However, their case study included large poloidal mag- 
netic fields, B > 10 12 G in both core and crust components, 
which would tend to violate the observed spin-down limit of 
B s < 2.0 x 10 11 G for PSR J0821-4300. 

5 . COMPARISON WITH OTHER PULSARS 

The ultimate goal of this field is to infer the equation of 
state and measure the radius of the NS. Some progress on 
these fronts has been made with high-quality data from mil- 
lisecond pulsars (MSPs). Using an unmagnetized hydrogen 
atmosphere model fitted to the spec tra and pulse profile o f the 
nearest known MSP J0437-4715. lBogdanov et al.l (120071) de- 
rived 6.8 < R < 13.8 km (for M = 1AM®). IBogdanov etall 
(2008, 2009) also obtained lower limits on R modeling X-ray 
observations of MSPs J2124-3358 and J0030+0451. Black- 
body emission was not abl e to fit the pulse profile s, thus re- 
quiring a NS atmosphere. Bogda nov et al.l d2007l) assumed 
an identical pair of polar caps, but fitted two temperatures to 
each, as required by the data, which can be understood as non- 
uniform heating by backflowing particles from the magneto- 
sphere givin g a concentric temperatu r e gradient, as originally 
modeled by iZavlin & Pavlovl (119981) . IBogdanov et"ai1 (12007b 
concluded that the magnetic dipole is not centered on the star, 
but must be offset by ~ 1 km to account for an asymmetry in 
the observed pulse profile. The data on PSR J082 1-4300 in 
Puppis A are not yet of a quality to search for such effects. On 
the other hand, the geometrical angles (£, ip) are not nearly as 
well constrained in the MSPs. (In the case of the binary MSP 
J0437-4715. it could be assumed that ip = 42° because that is 
the inclination angle of its binary orbit.) 

Pulsed light curves of the middle-aged pulsars PSR 
B0656+14, B1055-52, and Geminga, whose X-ray 
spectra are dominated by s urface thermal emission , 
have been model e d by iPagd (Q995), Page etliLl d!995l) . 
iPage & Sarmientol dl996l) . and lPerna etal.1 d2001l) . Begin- 
ning with RO S AT da ta it appeared that P SR B0 656+14 
dPossenti et alJ 1 1 996t iGreiveldinger et al.l 119961) . PSR 
B 1055-52 dOgleman & Finlevl I19931T and Geminga 
( Halpern & Ruderman 1993) had two thermal compo- 
nents with pulse-phase shifts of between 0.1 and 0.3 cycles, 
the hotter component being attributed to a heated polar cap. 
Follow-up observations at higher energy with ASCA found 
that th e harder components from Gemin ga ( Halp ern" & Wangl 
119971) and PSR B 1055-52 dWang et al.ll 1998b are better fitted 
by non-thermal power laws. As beamed emission from the 
magnetosphere, their hard X-ray pulses need not bear a sim- 
ple phase relationship to the soft thermal components. Only 
PSR B0656+14 continued to hav e two clear tempe ratures 
when observed at higher energy (Pav lov et al.ll2002ah . with 
only a weak nonthermal tail. 

Detailed study of the energy-dependent pulse profiles 
of these primarily thermal pulsars with XMM-Newton 
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dDeLuca et al.1 120051) confirm that PSR B0656+14 has two 
thermal components, with the hotter one interpreted as a small 
polar cap, shifted in phase by ~ 0.2-0.3 cycles from the softer 
emission. The spectrum of PSR B 1055-52 was fitted with 
two temperatures and a non-thermal power law, although it 
is difficult to explain why the hotter blackbody component 
has a pulsed amplitude of ~ 100 %. A case fo r a ho t po- 
lar cap on Geminga was made by ICaraveo et alj d2004l) and 
iDe Luca etaT] d2005l) , but lJackson & Halpernl d2005l) . analvz- 
ing the same data, did not find it to be necessary. Such a 
component does not make a significant contribution to the 
spectrum of Geminga at any energy, and its fitted pulse pro- 
file appears to have the same phase and similar shape as the 
power-law component, suggesting that it is a distinction with- 
out a difference. 

Despite these difficulties, the pulsed amplitudes and fitted 
areas of thermal X-ray emission from cooling neutron stars 
indicate that most have highly nonuniform surface temper- 
atures that may be regulated by their crustal magnetic field 
geometry. For example, the X MM-Newton observation of the 
middle-aged pulsar B 1706-44 dMcGowan et alj [2003) shows 
an asymmetric, double-peaked pulse profile whose T°° = 
8 x 10 5 K spectrum is compatible with the full NS area, while 
having a 22% pulsed fraction. In contrast, PSR J0538+2817 
appears to have only a hot polar cap of T°° = 2.2 x 10 6 K 
(McGowan et al. 2004). One of the most unusual results is 
the apparently thermal (T°° = 2.4 x 10 6 K) spectrum of the 
highfi -field (4.1 x 10 13 G) PSR Jl 1 19-6127 (tGonzdezetalJ 
2005), which has pulsations of amplitude 74% ± 14% that are 
only detected below 2 keV. 

Another family of thermally emitting p ulsars a re the 
nearby, isolated neutron s tars (INSs) dHaberll 120071; 
iKaplan & van Kerkwiikl 12009) with periods of 3.4- 11.4 s, 
and pulsed fr actions that range f rom 1.2% for RX 
J 1856.6-3754 dTiengo & Mereghettil [2007b to 52% for 
RX J1308.6+2127 dSchwope et al.l 120051) ~The latter authors 
fitted the double-peaked pulse profile of RX J 1308. 6+2 127 to 
a model of two small spots with temperatures of kTf° = 92 eV 
and fcr 2 °° = 84 eV separated by ~ 160° in phase. Timing 
measurements have revealed that these INSs have somewhat 
larger dipole m agnetic fields than most youn g pulsars, with 
B s > 10 13 G dKaplan & van Kerkwijkl 120091) . and may be 
significantly he ated by continuing magnetic field decay 
dPons et al.l 120091) . In this sense, they may be the ~ 10 6 year 
old descendants of magnetars. Several of their spectra have 
very broad absorption features that have been interpreted as 
ion cyclotron lines or, in the case of multiple features, pos- 
sibly atomic lines (Haberl 2007; Schwope et al. 2007). One 
of the best studied objects of this class, RX J0720.4-3125, 
shows a pulse phase shift of ~ 0.1 between soft and hard 
X-rays (Cropp er et al.l [2001), which suggests that there 
could be two spots of different temperatures. However, this 
interpretation is complicated by long-term (ye ars) changes in 
the shape of the spectru m and pulse profile (Ide Vries et al.l 
120041: iHohle et al] 12009). which lends support to the idea 
that, similar to the case of magnetars, localized and variable 
heating by magnetic field decay is responsible for relatively 



short-lived surface thermal structure. This is evidently the 
case for the transient magnetar XTE J1810-197, whose 
declining hot spot t e mpera t ures a n d areas were modeled by 
i Gotthelf & Halpernl d2005l 120071). iBernardini et al.1 d2009l) . 
Ulbano et al.1 d2009)) . and by iPerna & Gotthelfl d2008l) . who 
u sed a similar treatment as tha t presented herein. 

iPerez-Azorm et al.1 (l2006bh . IZane & Turolial (120061) and 
IZand d2007l) investigated models for INSs involving a combi- 
nation of star-centered dipole and quadrupole magnetic field 
components to explain their asymmetric pulse profiles. The 
properties of PSR J082 1-4300 in Puppis A may ultimately be 
ascribed to these same effects, but in a simpler system that is 
not variable in time. 

6. CONCLUSIONS 

We modeled the XMM-Newton spectra and pulsed light 
curve of PSR J082 1-4300 in Puppis A, one of three CCO 
pulsars whose dipole magnetic field strengths are measured 
to be less than those of all spin-powered pulsars of similar 
age (B s < 2.0 x 10 11 G in the case of PSR J0821-4300). The 
sizes and configurations of the surface hot and warm spots on 
PSR J0821-4300 are particularly well constrained. The two 
emitting areas differ by a factor of 2 in temperature and 20 in 
area, which conveniently endows them with similar luminosi- 
ties that fall in the XMM-Newton bandpass. The 180° phase 
reversal between the soft and hard X-ray pulse profiles reveals 
the antipodal geometry. It is especially significant that the X- 
ray spectra and pulse profiles of CCOs indicate considerably 
nonuniform surface temperatures. Many of the mechanisms 
that are held responsible for such effects in other classes of 
NSs are not expected to be operating in these anti-magnetars, 
which appear to be simple cooling neutron stars whose con- 
duction of heat from the interior is highly anisotropic. The 
essential problem in understanding CCOs is to explain how 
this is accomplished without creating a strong external dipole 
magnetic field. Our tentative hypothesis is that even CCOs 
have strong tangential fields buried in the crust that channel 
heat toward the magnetic poles, or external quadupole fields. 
But what is the geometry of that magnetic field? Although the 
orientation of the hot spots in PSR J082 1-4300 is determined 
to within 2°, the degeneracy of the model does not allow us 
to decide if the axis of the hot spots is nearly aligned with the 
NS spin axis or nearly perpendicular to it. Geometrical prob- 
ability, as well as the observation of larger pulsed fractions in 
other NSs, would suggest the latter. The actual geometry is 
probably fixed during the genesis of magnetic fields in these 
young NSs, which, in the case of CCOs, have not spun down 
since their birth and are likely to have preserved the natal B- 
field configuration. 

We thank the referee, Silvia Zane, for helpful comments 
on the paper. This work is based on observations obtained 
with XMM-Newton, an ES A science mission with instruments 
and contributions directly funded by ESA Member States 
and NASA. This work was supported by NASA XMM grant 
NN08AX71G. 
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